Chiral non-linear a-models as models for topological superconductivity 
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We study the mechanism of topological superconductivity in a hierarchical chain of chiral non- 
linear a-models (models of current algebra) in one, two, and three spatial dimensions. The models 
have roots in the ID Peierls-Frohlich model and illustrate how the ID Frohlich's ideal conductivity 
extends to a genuine superconductivity in dimensions higher than one. The mechanism is based on 
the fact that a point-like topological soliton carries an electric charge. We discuss a flux quantization 
mechanism and show that it is essentially a generalization of the persistent current phenomenon, 
known in quantum wires. We also discuss why the superconducting state is stable in the presence 
of a weak disorder. 



Introduction. In this paper we discuss models of topo- 
logical superconductivity, where fermionic currents are 
dynamically coupled with topological textures of chiral 
bosonic fields Q. The superconducting state in these 
models emerges due to an unbroken chiral symmetry. 
Electronic off-diagonal correlations in topological super- 
conductors are different from those of conventional BCS- 
type superconductors, but rather close to those of Lut- 
tinger liquids || . (i) Topological superconductivity takes 
place in standard models of current algebras (HH) (also 
called chiral non- linear tr-models), and (ii) the hydrody- 
namic action (and wave functions of low energy states) 
possesses a geometric phase (or $-term) - an imaginary 
part of the Euclidian action 



iSr = 2ir\, 



dx dt j 



N f 0. 



(1) 



Here j M = (vop,j) is an electronic current and density 
deviation from the equilibrium (we set vo — 1 in the 
following), Ai is the London's penetration depth, and 
Nf is a degeneracy of an electronic state (e.g., spin). 

The first term in ([!]) gives the London equations. ■& is 
an integer valued topological invariant. The theta term of 
(|l|) does not affect the equations of motion but the mon- 
odromy of the wave function. It is the geometric phase 
which carries an information about spacetime transfor- 
mation properties of excitations. In particular, it changes 
by 1 in the processes of a 27r-rotation of a group of Nf 
electrons around some axis or an interchange of positions 
of two groups of Nf electrons. In a multi-connected sam- 
ple, it also changes by 1 under the translation of Nf elec- 
trons along a noncontractable path. We show that this 
term provides a mechanism for flux quantization in a 
multi-connected geometry. The geometric phase is, per- 
haps, the most important difference from the s-wave BCS 
theory, where the Euclidian action is real. It emphasizes 
a quantum interference aspect of the mechanism. While 
discussing the z9-term, we rely on the results of Ref. || . 

Models of current algebras. The models of topological 
superconductivity in two and three spatial dimensions 



are the standard non-linear cr-models of current algebras 
(see e.g. [pi). These models describe Dirac fermions in- 
teracting with a chiral boson field 



Co = 



iD - mn ■ f J ijj + W 2 [n], (2) 
iD - m(ir + i-f 5 jr ■ f) ty + W^n], (3) 



where iD = ^^(id^ + A^) and 7 M are 2x2 and 4x4 
Dirac matrices in spatial dimensions two and three. The 
fermions are SU(2) doublets and f = (t\, to,, T3) are the 
Pauli matrices. In addition, a fermion variable tp carries 
a (suppressed) flavor index running from 1 to Nf. The 
parity of the number of flavors Nf plays an important 
role in what follows. The chiral fields take values on two 
and three-dimensional spheres thus matching the number 
of spatial dimensions: ft 2 — 1 and 7r 2 + fr 2 = 1. We test 
the system by an external electromagnetic potential A^. 
Finally, W2. 3 describe a bare action of the chiral field, 
other than the dynamics induced by fermions. In 2D, 

e.g., it is W 2 = H (ft - s 2 (djft) 2 ) +■■-. The models have 
an electronic origin and can be viewed as an effective 
description of certain interacting electronic systems. We 
do not discuss that here. 

We emphasize that the essential properties of these 
models are inherited from the model of Frohlich's ideal 
conductivity in one spatial dimension j|. We, there- 
fore, add the latter to the hierarchy. In this model, ID 
fermions interact with a chiral field taking values on a 
circle: 



iD 
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It is well known 0] that 27r-phase solitons of the Frohlich 
model are charged and move freely through the system, 
making it an ideal conductor. It is not a superconductor, 
however, since the restricted spatial geometry does not 
allow for a transverse motion of a charge flow. As a result, 
a weak disorder or commensurability effects destroy the 
ideal conductivity |11|. We will see that in dimensions 
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greater than one (models (|2|||)) a phenomenon similar to 
Frohlich's ideal conductivity results in superconductivity 
even in the presence of weak disorder. 

The physics of the models is controlled by the 

chiral symmetry. In ID it is: 4> — * <fi + a, ijj — > e~^ lba ip. 
In 2D it is: nr — > U^n-fU and ip — > U -1 ^, where U is an 
SU(2) matrix. In 3D it is g -> U~ x gXJ and t/» -» C/" 75 ^, 
where g = 7ro + ittt. 

We argue that if the global chiral symmetry is unbro- 
ken, the ground states of the 2D and 3D models are su- 
perconducting. These models then exhibit a Meissner 
effect and a flux quantization with period if Nf is 

even and if Nf is odd, where $o = hc/e is the flux 
quantum. Notice that the period is always a fraction of 
$o with an even denominator regardless of whether the 
number of fermionic species is even or odd. For the most 
interesting cases of polarized (Nf — 1) and unpolarizcd 
(Nf = 2) electrons, the flux is quantized with the period 
$ /2. 

Electronic and topological currents. In every model (||- 
||) a coordinate space matches the target space of chiral 
fields. Therefore, there exist spatial configurations of a 
chiral field which wrap over the target space with an 
arbitrary integer winding number Q - topological charge. 
They are solitons: phase solitons in ID and skyrmions in 
D = 2,3. The topological charge Q may be written as 
a spatial integral Q — J d D x Jq of the zeroth component 
of topological current J M : 

= 7r£tiud v <i>, (5) 

ZTT 

J n = -^-tuuxn ■ d v n x d\n, (6) 
on 

Jfl = 127T 2 € ^\piabcdn a d v TTbd\TT c dpTTd- (7) 

In configurations without singularities, a topological cur- 
rent conserves d^J^ = 0. This reflects the fact that topo- 
logical charge cannot be changed by smooth variations of 
chiral fields. 

It is known that the energy spectra of the models (|[]^) 
in a smooth static soliton background with a topological 
charge Q have an additional NfQ levels with negative en- 
ergy, compared to the spectrum in a topologically trivial 
background. Being filled, these levels create an electric 
charge localized on the solitons. The motion of the soli- 
tons then results in an electric current since the electric 
charge moves together with a soliton as long as the mo- 
tion is adiabatic [l2| . This is summarized by 

fj d D x = N f Q; j fl (x)=NfJ^(x) + ---, (8) 

where j„ = (V>7 M '0} is an electronic current in a back- 
ground chiral field and J M is a topological current of a 
chiral field f^fy- The density of the electronic current 
(H) has corrections (suppressed by the gap m) due to 



terms with higher gradients. However, the corrections 
do not have a topological nature. The change of total 
charge in a soliton background is exact. 

Spin and momentum of a soliton and the fl-term. Intu- 
itively, it is appealing that a soliton carrying a fermionic 
charge NfQ has a spin NfQ/2 and statistics that depend 
on whether NfQ is even or odd. It has been shown in 
Refs. Jl3|,^] that this is indeed true - an intrinsic quan- 
tum number of a soliton - a fermionic number (||) con- 
verts into rotational and translational quantum numbers. 
Spin and statistics of solitons are described by a geo- 
metric phase of the wave function. The geometric phase 
changes by nNfQ under an adiabatic 27r-rotation of a 
soliton or an adiabatic exchange of the positions of two 
solitons. In a multi-connected geometry, the geometric 
phase also changes by the same amount under an adi- 
abatic translation of a soliton along a noncontractable 
cycle. Along such a path, the wave function acquires a 
phase e l7rJV/< 2, which indicates the spin and statistics of 
a soliton. The geometric phase also results in a shift in 
quantization of the orbital momentum L = Nf ® + I and 
momentum IT = j^(Nf® + ni) along a noncontractable 
cycle (here I and are integers and Li is the length of 
the cycle). 

We write fermionic parts of our models as C — ijjUip 
with T> an operator, dependent on chiral fields. Then the 
induced action W of the non-linear er-model is a result of 
an integration over fermions in a fixed background con- 
figuration of a chiral field e~ w — J e~ § cdxdt dipd4> — 
(DetD) f. Then the geometric phase appears as an 
imaginary part of the Euclidian action of a non-linear a- 
model. Indeed, in Refs. jl3],[§ it has been shown that the 
geometric phase in the models (||-||) is equal to NfW times 
a topological invariant $ that characterizes a topology of 
spacetime configurations of the chiral fields <p, n, (ttq,tt). 
The $-term is also referred to as a global chiral anomaly, 
distinct from the local anomaly given by (||). 

The $-term depends on the compactification of the Eu- 
clidian spacetime. If it is a sphere, the ^-invariant is a 
representation of the homotopy groups ttd + i(S d ). It is 
zero in ID, any integer (the Hopf number) in 2D, and 
or 1 in 3D. 

In a more physical case, the Euclidian time is a circle 
and the coordinate space may also have some noncon- 
tractable cycles. The topological i9-term then is corrected 
by terms describing the motion of solitons along the cy- 
cles. The explicit expression for the z9-term requires a 
higher profile of the topology |Q , but its physical mean- 
ing is simple. E.g., an adiabatic rotation of a soliton by 
an angle (j> results in the t9-term = Q J ifij^ (see || 
and references therein) . If a soliton is adiabatically trans- 
lated along a noncontractable spatial cycle, the $-term is 
fl[4>] — Q J Ut^-, where u is a displacement of the soliton. 
The latter is the only term which remains in ID, where it 
is known to be responsible for the parity effect in the per- 
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sistent current phenomenon observed in quantum wires 

Non-linear a-models. Fermionic number (H) and the 
i?-term determine the imaginary part of the Euclidian 
action of a non-linear cr-model. Its real part is given by 
a regular gradient expansion. To leading order in 1/m, 
we have (see e.g., Q and references therein) 

^ = inti[<f>] + J d 2 x (iA^ + ^(d^ , (9) 
^ = i*d[n] + J d 3 x (iA^ + £ (^n) 2 ) , (10) 
^ = M\g] + J d 4 x (iA^ + ^tr^S -1 ^)) . (H) 

where F 2 = -J-^rn 2 In — with A being an ultraviolet cut- 
off and g = ttq + inr a unitary matrix. The first two 
terms represent global and local chiral anomalies. 

London's hydrodynamics. The hydrodynamics of 
a superconductor is defined by a rigidity with re- 
spect to transverse charge currents (j[i(k), j v (—k)) = 
(47rA 2 e 2 )~ 1 (8 l _ ll , — ty t ;" ) at k — > and can be summa- 
rized by the London Lagrangian: C = 2nXj j j 2 . 

In topological liquids, charge current is equal (at small 
k) to the topological current (||), so that a criterion for 
superconductivity translates to a transverse rigidity of 
topological currents 

(J,(fc), J v (-k)) = (4?rN 2 \ 2 L e 2 )-i (s^ ^) . (12) 

In ID the first terms of the gradient expansion (^ al- 
ready give ( |l2"| ) . We will argue that the transverse rigidity 
of topological currents takes place in any dimension if the 
global chiral symmetry is unbroken, i.e., if the chiral fields 
are in a quantum disordered state: (ft(x) , ft(0)) — * at 
large x in 2D, or (g(x), <? _1 (0)) — ► in 3D (it is always 
true in ID model). 

The situation in higher dimensions is obscured by the 
fact that the models (||||) ar e not renormalizable and 
that the symmetry of the ground state depends on the 
effect of small distances not included in the continuum 
effective model. To implement the state with unbro- 
ken chiral symmetry (a quantum disordered state), one 
might choose a non-electronic part W2.3 of the models 
(HH). The latter does not contribute to the anomalous 
terms. In addition, we assume that the correlation radius 
of the quantum disordered state r c is large so that the 
chiral field is smooth and the electronic density follows 
the topological density according to eq. (||) . 

The arguments that a quantum disordered state of a 
non-linear cr-model, regardless of its dimension, is rigid 
with respect to topological excitations, are general |16|. 
Here we adopt them for the 2D model following Ref. [p|. 
We first write the non-linear cr-model (JIG) in the CP 1 
representation by writing (d^ri) 2 as Az*{— id^ — ap) 2 z, 



where z = [z\,z%) is a complex doublet with constraint 
z*z = 1 related to n as ft — zWz. The field strength 
of the gauge field is proportional to the topological 
current of ft according to J M = ^e^ vX d v a\ - it describes 
the fluctuations of the solid angle formed by the n-field. 
The rigidity of topological currents (|l^) means that the 
propagator of the gauge field a M induced by z-bosons is 
massless. 

It is obvious that if the chiral field n is ordered, than z- 
bosons are condensed and ( J^(fc), J v {— k)) ~ nk 2 , where 
the stiffness k is proportional to the amount of the con- 
densate of z-bosons. This combined with (|^) leads to an 
insulating behavior. 

On the contrary if the chiral symmetry is unbroken 
the system does not respond to the constant transverse 
field a M . This means that the gauge field is massless 
(a u (k), a v {— k)) ~ (r c fc 2 ) -1 , which leads to ( |l2| ) with 
r c ~\le 2 0k 

Combining (|12|), the fermionic number of solitons, and 
the geometric phase (|^-[ll]), we obtain the hydrodynamic 
description of a topological superconductor ([j]) in D=2,3 
and Frohlich's ideal conductor in D=l. In terms of topo- 
logical currents, it is 

S = N f J {2nN f X 2 L e 2 J 2 + ieJ^) dxdt + wriV/tf. (13) 

Disorder. It is generally correct that if the Meissner 
effect takes place in a system with no disorder, then a 
superconducting state remains intact in the presence of 
weak disorder, regardless of the mechanism that led to 
superconductivity. In ID there are no transverse direc- 
tions and no true Meissner effect. As a result the Frohlich 
ideal conductivity is destroyed by the backscattering of a 
single impurity [ pT| . A 2fc/-harmonic of the impurity po- 
tential Vzkt adds the term \V2kt I cos(c6— a) to the effective 
Lagrangian (^), where a is the phase of the potential. It 
pins a soft translational mode of the charge density wave 
at c6 = a, makes fluctuations of <f> massive, and suppresses 
correlations of topological currents at small momenta. 

In higher dimensions, the existence of transverse cur- 
rents eviscerates this effect. In 2D, similar to the ID 
case, we add a backscattering term V^i-.n to (|Io[), where 

Vski is a vector, representing a 2k f harmonic (scattering 
between Dirac points of the spectrum) of an impurity 
potential. This term has a tendency to align the n-field 
along Vik s , but does not change the correlation of topo- 
logical currents (solid angles of ft) as long as the potential 
is smaller than a scale provided by the correlation radius 
r c . The system remains in the quantum disordered phase, 
while the superfluid density is reduced by the amount of 
the order V 2 kf . One can see this formally by writing 

the impurity term in CP 1 representation V^kj^fz. This 
term shifts the local chemical potentials of the z bosons in 
opposite directions for the bosons with isospin along and 
opposite to V2k f ■ If z-bosons are far from the condensa- 
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tion point (ordering of ft), the small change of chemical 
potential cannot drive the condensation, and fluctuations 
of an internal gauge field remain massless. The super- 
conducting state cannot be destroyed by an infinitesimal 
amount of impurities. The important difference with ID 
is that the impurity potential does not couple to a topo- 
logical density of solitons (i.e. to solid angle of n-field in 
this case) . Similar arguments are valid for the 3D model. 

Flux quantization. The geometric phase mNf& in the 
low energy action (Q) provides the mechanism of the flux 
quantization in topological superconductors in a multi- 
connected geometry. Consider a sample in the shape of a 
torus, threaded by magnetic flux <f> through some partic- 
ular torus cycle, which we label by 1. Due to the gauge 
invariance all energy levels are <E>o-periodic functions of 
$. We will show that the energy of a macroscopic sample 
at a fixed chemical potential is periodic with a smaller 
period: $a/Nf if Nf is even and <&o/2 Nf if Nf is odd. 

Starting from the low energy action fl3|), we pass to 
the Hamiltonian in a sector with a given topological 
charge Q. We introduce displacements u according to 
Jo = -Vu, J = u. The relevant part of the $-term 
can be written in terms of zero modes (the homogeneous 
part) of displacements M, = J UidP x. E.g, a soliton mov- 
ing around the cycle 1 contributes •& = Q J j^dt (similar 
to the one in ID, where u\jL\ = c/)/2tt). 

Neglecting all excitations transverse to the torus cy- 
cle pierced by the flux, we set all currents except for 
Ji to zero. Then we neglect all harmonics of the dis- 
placement u\ except for the zero mode. Their contri- 
bution to the energy does not depend on the flux. Set- 
ting the gauge potential to be a constant along the cycle 
1, A = (27t$/Li$ 0j 0, 0), we obtain the flux dependent 
part of the Lagrangian (we return to the real time) 

Varying it with respect to u\ we find canonical momen- 
tum Ei = 4^V 2 A|e 2 ^ - t^(§ + £), where we ob- 
serve a shift of the momentum by the z9-term. Repre- 
senting the canonical momentum of the zero mode as 
III = — i-S—, we obtain the wave function of the zero 
mode exp (z 2 ^" 1 Ux), where n\ is an integer, and the flux 
dependent part of the energy of the sliding charge density 
mode 

If Nf is even, it is obvious that the energy is a $>o/Nf- 
periodic function of the flux (c.f., ||). However, the en- 
ergy remains $o/2-periodic even at Nf = 1. In fact, if 
Nt is odd the energy is <I>o/2iVf -periodic. In macroscopic 
samples, the energy of vacuum polarization by a soliton 
- the last term in (jTJ) - is negligible compared to the 



kinetic energy. If Nf is odd, the ground state possesses 
no solitons Q = as long as <f> < <&a/ANf. If, however, 
<f>o/4iV/ < $ < $ /2Nf, the ground state at the fixed 
chemical potential acquires a topological charge Q = 1 
and Nf additional particles. The freedom to create a 
soliton, when the flux exceeds §>o/4Nf, halves the naive 
period $>o/Nf of the kinetic energy. 

In the ID setup (a quantum wire) V = L\, \\e 2 = 
l/4Nf, and ( |l4| ) gives the persistent current j — 
—dE/d& for the Peierls-Frohlich model (Q). In this case 
the vacuum polarization energy caused by a static soliton 
(the last term in eq.(p^)) destroys periodicity for odd Nf 
(c.f., H). 

This mechanism is to be compared with the familiar 
flux quantization in the BCS theory. In that flux- 
dependent part of the energy is given by zero modes of 
the phase (p — 2im-^ of the order parameter along a 

noncontractable cycle: E n — 2 \^ 2 l 2 (f + 5^) " ^ ^ s 
$o/2-periodic and has deep minima at fluxes which are 
multiples of $0/2. Although the topological mechanism 
leads to a similar result, the state counting and the un- 
derlying physical picture is different. 
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